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We study the link overlap between two replicas of an Ising spin glass in three dimensions using the 
Migdal-Kadanoff approximation and scaling arguments based on the droplet picture. For moderate 
system sizes, the distribution of the link overlap shows the asymmetric shape and large sample- 
to-sample variations found in Monte Carlo simulations and usually attributed to replica symmetry 
breaking. However, the scaling of the width of the distribution, and the link overlap in the presence 
of a weak coupling between the two replicas are in agreement with the droplet picture. We also 
discuss why it is impossible to see the asymptotic droplet-like behaviour for moderate system sizes 
and temperatures not too far below the critical temperature. 

PACS numbers: 75.50.Lk Spin glasses 



I. INTRODUCTION 



A series of computer simulations performed in the past 
few years [|l]J^] appear to support the claim that the 
three-dimensional Edwards- Anderson spin glass shows 
signatures of replica-symmetry breaking (RSB), imply- 
ing the existence of infinitely many pure states in the 
low-tempcrature phase. In contrast, almost rigorous ar- 
guments 1 3 ] and recent experiments Q favor the droplet 
picture [p - 7[ with only one pair of pure states. 

In a recent paper Q , we have suggested that due to the 
high temperatures and small system sizes the computer 
simulations are strongly affected by the critical point 
and do not reflect the true low-temperature behaviour. 
This suggestion was supported by a numerical calculation 
of the Parisi overlap function using the Migdal-Kadanoff 
approximation (MKA) . For the system sizes and temper- 
atures typically used in computer simulations we found 
overlap functions similar to those in [Q , however for lower 
temperatures we found agreement with the predictions of 
the droplet picture. That the results of computer simu- 
lations are strongly affected by the critical point can also 
be concluded from Q , where the Parisi overlap function 
shows critical scaling (with effective exponents) down to 
temperatures ~ 0.8T C . 

In recent publications Jl0|-4l2|, it is claimed that non- 
trivial behaviour of a quantity called the link overlap is a 
reliable indicator of RSB. However, in order to place such 
a claim on solid ground, one would have to show that the 
data cannot be interpreted within the framework of the 
droplet picture. It is the goal of this paper to furnish this 
discussion which has been missing so far. As in || , we use 
the MKA which is known to agree with the droplet pic- 
ture. Our results show, as in the case of the Parisi over- 
lap function, that several nontrivial features attributed 
to RSB are in fact due to finite-size effects, and that the 
numerical data on the link overlap published so far are 
indeed in agreement with the droplet picture. We also 
derive expressions for the effective coupling at any tem- 
perature as a function of system size and find that one 



indeed needs rather large systems or low temperatures to 
see droplet-like behaviour. 

The outline of this paper is as follows: After introduc- 
ing the model and defining the quantities to be evaluated, 
we present first our analytical and numerical results for 
the link overlap distribution function. Then, we evalu- 
ate the link overlap in the presence of a weak coupling 
between the two replicas. In the following section we 
explain why finite size effects are so large for the three- 
dimensional Ising spin glass. Finally, we summarize and 
discuss our findings. 



II. DEFINITIONS 

The Hamiltonian Ho of the Edwards-Anderson (EA) 
spin glass in the absence of an external magnetic field is 
given by 



@H = - J v a,o 

(id) 



j > 



where (3 = 1/ksT. The Ising spins can take the values 
±1, and the nearest-neighbour couplings Jy are inde- 
pendent from each other and gaussian distributed with a 
standard deviation J. 

It has proven useful to consider two identical copies 
(replicas) of the system, and to measure overlaps between 
them. This gives information about the structure of the 
low-temperature phase, in particular about the number 
of pure states. The quantity considered in this paper is 
the link overlap 



(e) = (l/N L )J2(WjnTj) 



(1) 



where the sum is over all nearest-neighbour pairs of 
a lattice with Nr, bonds and N sites, and the brackets 
denote the thermal and disorder average, a and r denote 
the spins in the two replicas. The Hamiltonian used for 
the evaluation of the thermodynamic average is 
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PH[<t,t]=PH o [j]+0H o [t] 



e^ViCTjTiTj 



(id) 



(2) III. THE PROBABILITY DISTRIBUTION OF 

THE LINK OVERLAP 



where Hq is the ordinary spin glass Hamiltonian given 
above, and the term in e introduces a coupling between 
the two replicas. 

In cases where the random couplings are taken to 
have the values ±1, the link overlap is identical to the en- 
ergy overlap. The main qualitative differences between 
the Parisi overlap 
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and the link overlap are (i) that flipping all spins in one 
of the two replicas changes the sign of but leaves 
g( L ) invariant, and (ii) that flipping a droplet of finite 
size in one of the two replicas changes q( p ' by an amount 
proportional to the volume of the droplet, and q( L * by an 
amount proportional to the surface of the droplet. 

Below, we will show that, just as for the Parisi overlap, 
the MKA can reproduce all the essential features of the 
link overlap found in Monte Carlo simulations. These 
results refute the claim made in |12| that the agreement 
between the MKA and simulations for the Parisi overlap 
reported in || is a mere coincidence that does not extend 
to the link overlap. The conclusion must be drawn that 
there is no evidence for RSB in three dimensional Ising 
spin glasses. 

Evaluating a thermodynamic quantity in MKA in 
three dimensions is equivalent to evaluating it on a hierar- 
chical lattice that is constructed iteratively by replacing 
each bond by eight bonds, as indicated in Fig. The to- 
tal number of bonds after / iterations is 8 7 , which is iden- 
tical to the number of lattice sites of a three-dimensional 
lattice of size L = 2 1 . Thermodynamic quantities are 
then evaluated iteratively by tracing over the spins on 
the highest level of the hierarchy, until the lowest level 
is reached and the trace over the remaining two spins is 
calculated [Q. This procedure generates new effective 
couplings, which have to be included in the recursion 
relations. In pj , it was proved that in the limit of in- 
finitely many dimensions (and in an expansion away from 
infinite dimensions) the MKA reproduces the results of 
the droplet picture. We have shown in Q that the MKA 
agrees with the droplet picture in three dimensions as 
well. For this reason, no feature that is seen in MKA can 
be attributed to RSB. 



FIG. 1. Construction of a hierarchical lattice. 



We first set the coupling strength e in Eq. (|2|) to zero 
and study the probability distribution P(q^) of the link 
overlap, averaged over a sufficiently large number of sam- 
ples. RSB should manifest itself in P(q^) according 
to ]Trj| , pT| in an asymmetric (non-Gaussian) shape and a 
nonzero width even at infinitely large system sizes. Fur- 
thermore, the link overlap for single samples should show 
large variations between different samples. Here, we show 
that the asymmetric shape and large sample to sample 
variations can even be seen in MKA for moderate sys- 
tem sizes and can therefore not be taken as evidence for 
RSB. The only reliable indicator for RSB would be a 
width of P(q( L >) that does not shrink with increasing 
system size. However, the only Monte Carlo simulation 
data published so far for P(q( L ') [[llj are taken for a four- 
dimensional Ising spin glass in a magnetic field, and they 
show a shrinking width, in agreement with the expecta- 
tions from the droplet picture. 

We have obtained the function P{q^) by first calcu- 
lating its Fourier transform, 



F{y) = ( exp liy^2 ■ 



(v) 



The coefficients a n in 



N L /2 

P(q {L) )= J2 a n 5(q^ ~2n/N L ). 

n=-N L /2 

are then found from F(y): 

rnN L /2 

a n = (l/nN L ) F(y)exp(2iyn/N L )dy. 

J-TrN L /2 

Figures^, ^| and |2| show our result for P(q( L >) in MKA 
for three different temperatures. All curves have been 
averaged over several thousand samples. 




FIG. 2. The link overlap distribution P(q {L) ) for T 
and L — 4,8,16,32 (from widest to narrowest curve). 
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FIG. 3. The link overlap distribution P{q (L) ) for T = 0.7T C 
and L — 4,8,16,32 (from widest to narrowest curve). 
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FIG. 5. The link overlap distribution P(</ L ') for four dif- 
ferent samples at L — 8 and T = Q.7T C . 
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FIG. 4. The link overlap distribution P(q (L) ) for 
T = 0.45T C and L — 4,8,16,32 (from widest to narrowest 
curve) . 
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FIG. 6. The link overlap distribution P(q ( - L ^) for three dif- 
ferent samples at L — 16 and T = 0.7T C . 



The curves for small L are asymmetric at T = T c with 
a tail on the right-hand side. With decreasing tempera- 
ture, the asymmetry becomes stronger, the tail moves to 
the left-hand side, and a shoulder is formed. All these 
features seem to be finite-size effects, as they become 
weaker with increasing system size. 



In particular for L — 8, there are large variations be- 
tween different samples, a feature that is usually assumed 
to be a clear indicator of RSB. However, since the MKA 
does not show RSB, we must asign this feature to finite- 
size effects. 

Finally, let us study the width 



(q( L ) - q( L r) 2 P(q( L r)dq( L ) 



Figures | and | show P(q^) for single samples at 
T = 0.7T C and for L = 8 (Fig. §) and L = 16 (Fig. g). 



ofP(q( L )). Fig. @ shows our results on a double logarith- 
mic plot, together with a power-law fit q( L *> ~ L~ w . 
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FIG. 7. Width of P(q <L) ), A<j (L) , as a function of L for 
T — T c (circles), T = 0.7T C (squares) and T = 0.45T C (dia- 
monds), together with a power-law fit. The exponents of the 
three fits are 1.43, 1.33 and 1.23. 

At T = T c , the exponent lu is sufficiently close to 
d/2 = 1.5 to suggest that the leading contribution to 
Ag( L ' comes from the superposition of independent con- 
tributions of the different parts of the system, just as it 
does for higher temperatures. Below, we shall see ex- 
plicitely that critical point nonanalyticities are indeed 
subleading to the regular contributions. 

Within the framework of the droplet picture, the value 
of to at the zero temperature fixed point can be calculated 
as follows: The main excitations at low temperatures or 
large scales are droplets of flipped spins of a radius r ~ L 
in one of the replicas. Such droplets occur with a proba- 
bility proportional to 

{L/r) d r- 6 kT, 
and each of them makes a contribution of the order 

r 2d„L-2d 

to (Aq( L )) 2 . d s is the fractal dimension of the droplet 
surface and is d — 1 in MKA, and 6 is the scaling ex- 
ponent of the coupling strength and has a value around 
0.24 in MKA in d = 3. We therefore find 

Aq^ ~ VkTL d *- d - e / 2 , 

giving uj ~ 1.12 in MKA. Our above data for T = 0.45T C 
are not far from this result. For larger system size, they 
must ultimately converge to it. Just as in the case of 
the Parisi overlap ||, the crossover from critical to low- 
temperature behaviour is so slow that even for T ~ 0.4T C 
the asymptotic regime is not reached for system sizes up 
to 32, and the curves appear to show effective exponents. 

For a cubic lattice, we have d s ~ 2.2 and 9 ~ 0.2, 
predicting a value lu ~ 0.9 within the framework of the 
droplet picture. This has to be compared to the RSB 
scenario, where u = 0. In Q, it is claimed that Monte- 
Carlo simulation data at T ~ 0.6T C and L < 12 show 



already the signatures of RSB. However, as yet there is 
no published data for Aq( L \ Only if a value of lu smaller 
than 0.9 is found, one can conclude that the droplet pic- 
ture is inappropriate and that RSB occurs. As long as 
lu appears to be larger than 0.9, the data are compatible 
with the droplet picture and are affected by finite-size 
effects. 



IV. THE LINK OVERLAP IN THE PRESENCE 
OF A COUPLING BETWEEN THE TWO 
REPLICAS 

In p0| , ^2[ , the authors suggested studying the expec- 
tation value of the link overlap in the presence of a cou- 
pling between the two replicas, q^ L \e), in order to test 
whether a system shows RSB. A positive coupling e in 
Eq. (|^) favours a configuration where both replicas are in 
the same state, while a negative e favours configurations 
with smaller overlaps. If the RSB scenario were correct, 
the distribution P(q^) would have a finite width even 
for L — > oo and range from some minimum value q m i n to 
a maximum value q ma ^. Consequently, the expectation 
value of the link overlap, q( L \e), would have a jump from 
<7min to <7 m ax at e = in the thermodynamic limit L — > oo 
]12| ■ In contrast, as we will show next, the droplet picture 
predicts a continuous and at e = nonanalytic function 
q^ L \e) [jl5|]. Part of these results were published in 

Within the droplet picture, the scaling dimension of e 
in the spin glass phase can easily be obtained. At length 
scale 1, e is equivalent to the energy cost of flipping one 
spin in one of the two replicas. On a scale I, this becomes 
the energy cost of flipping a droplet of radius I in one of 
the replicas, which is proportional to el ds . The scaling 
dimension of e is therefore d s , where d s is the fractal di- 
mension of the droplet surface. Equivalently, d s is the 
fractal dimension of a domain wall. Within MKA, we 
have d s = d — 1. The same value d — 1 for the scaling 
dimension of e is also obtained from an analytical calcu- 
lation of the recursion relations for e and the strength J 
of the random couplings near the T = fixed point. 

The positive dimension of e implies that the coupling 
between the two replicas is a relevant perturbation and 
that on large scales a behaviour different from that of an 
independent system can be seen. When e is positive, the 
energy cost el ds for the excitation of droplets of radius I 
leads to the suppression of droplets larger than 

I* - (kT/e) 1/d \ 

On scales beyond I*, droplet excitations must occur in 
both replicas simultaneously. This costs twice the en- 
ergy of a single droplet in an independent system. The 
coupled system thus behaves on large scales exactly like 
a single system with twice the coupling strength J. 

As stated in the preceding section, for e = 
droplets of size / occur with a probability proportional 
to (L/l) d l~ e kT in a system of size L. Since for posi- 
tive e droplets of size greater than I* are suppressed, the 
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change in the link overlap due to a small positive e can 
be written as 

L 

q( L \e)-q^(0)~ hT^ 1 ^'"- (3) 

i>l* 

In order to suppress each droplet only once, the sum 
must be taken over distinct length scales I, e.g., I — 
I* ,21* , and it is proportional to the first term. 
Therefore, we can write ( e ) ~ (°) - kT{l*) d ^ d - e , 
and using the expression for I* given above we have, 

q {L \t) - q {L) (0) ~ kT(e/kT)( d+e - d ° )/d \ (4) 

For negative e, flipping a droplet of radius / in one of 
the replicas changes the system's energy by an amount 
proportional to I s — |e|Z ds , which is negative for I > l c 
with l c ~ \e\ y /^- d '\ Therefore, there is a proliferation 
of droplets beyond this length scale and the spin glass 
state is completely restructured. We followed the flow of 
the parameters e, J, and Ae (the width of the distribu- 
tion of e) under a change of scale in the MKA and found 
that Ae diverges, while J and e eventually decrease to 
zero. Such a system is an Edwards- Anderson spin glass 
with the effective spins pi — a{Ti. 

Since droplets of size larger than l c proliferate for neg- 
ative e, the change in the link overlap is given in this case 
by 

jW(e) - g (i) (0) ~ -(l B ) d '- d ~ -| c |(«*-«*-)/C«'--e). ( 5 ) 

We thus find that q {L \e) - q^(0) has the form 
^±l e | ' with values A and A that depend on the sign of 
e. Within MKA, it is A+ ~ 0.62, and A_ ~ 0.57. For a 
cubic lattice, A+ ~ 0.45, and A_ ~ 0.40. 

For finite temperatures and small systems, there are 
corrections to this asymptotic behaviour due to finite-size 
effects which replace the nonanalyticity at e = with a 
linear behaviour for small |e|, and due to the influence 
of the critical fixed point, where the leading behaviour 
is linear in e (see below). As we have argued in S, the 
influence of the critical fixed point changes the apparent 
value of the low-temperature exponents for the system 
sizes studied in Monte-Carlo simulations and the MKA. 
The data shown in p2| with an apparent value of 0.5 
for A± arc fully compatible with the above predictions of 
the droplet picture. There is no indication of a jump at 
e = in gW(e), which would be the signature of RSB. 
For the MKA, the apparent exponent at 0.7T C is close to 
I for L ~ 16, leading to the "trivial" behaviour found in 
[ p2[ . However, at lower temperatures, for the same small 
system sizes the above-mentioned nontrivial features pre- 
dicted by the droplet picture become clearly visible, as 
shown in Fig. 0. 
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FIG. 8. q (L) (e) in MKA at T ~ 0.38T C (bottom) and 0.14T C 
(top) as function of sign(e)|e| ' 6 for various system sizes. Each 
curve is averaged over several 1000 samples. 

Let us conclude this section with a discussion of the 
link overlap at the critical temperature. Fig. ^| shows our 
result in MKA. 
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FIG. 9. q [L) {e) in MKA at T ~ T c for various system sizes. 

Clearly, the curves are linear at e = 0, indicating that 
the regular part dominates over the singular, critical con- 
tribution. This conclusion is confirmed by studying the 
scaling dimension of e at T c . Iterating the recursion re- 
lations for the coupling constants, we find that under a 
change in length scale, x — > x/b, we obtain e — > b^e with 
<f) ~ 1.14. Now, q [L) {e) can be obtained from the free 
energy via the relation 

q (L) (e) = (l/3N)d(lnZ)/de, 

implying a scaling behaviour q^ — > b d ~^q^ L \ Substi- 
tuting b with e gives then the relation 

q( L ) ^ e (<i-0)/0 ^ e l-6 

Compared to the linear regular part, this singular depen- 
dence cannot be seen for small e. 
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V. FINITE SIZE EFFECTS 

As we have argued throughout this paper and in earlier 
work H [^5| , finite size effects appear to be large for the 
Ising spin glass in three dimensions. To understand this 
behaviour, we iterated the MK recursion relations start- 
ing at various temperatures below T c . In Fig.jio] the solid 
lines show the data for J 2 as a function L = 2" where 
n is the number of iterations of the MK recursion rela- 
tions. For large L, one expects J 2 ~ 2 2n6 . However, it is 
apparent that even at T ~ 0.7T C one needs system sizes 
L ~ 100 to see this behaviour. Because the change in 
the slope din J 2 /d\nL is so slow, J 2 (L) appears to be de- 
scribed by a power law with some effective exponent over 
small windows of one decade in L, just as we found for 
P( p )(0) in ||. To understand the large crossover regime, 
we consider an expansion around the zero-temperature 
fixed point, where the effective temperature T — 1/J at 
a length scale L can be written as 



dT/d\nL = -9T + AT 3 + BT 5 



(6) 



where A, B, ... are constants and even order terms are 
absent because T — ► — T (or { Jy} — > {— Jij}) is a symme- 
try of the Hamiltonian. Now, d — 3 is close to the lower 
critical dimension for Ising spin glasses (so that T c is in 
some sense small) , and one might expect that truncating 
the above equation after the first few terms will give a 
good approximation for T or equivalently J throughout 
the low temperature phase. We now show that this is 
indeed the case: keeping terms up to T 5 gives a good 
description of the MK data of Fig. 10|. 

The analysis is as follows: At T c , dT/dL — so that 
9 = AT 2 + BT C A . For small deviations away from T c i.e. 
T = T c + 5T, the correlation length exponent v is defined 
via d\n(ST)/d\nL = l/v, leading to l/v = 26 + 2BT*. 
Solving for A and B in terms of v 1 9 and T c , we find 
A = (29 - l/2v)/T 2 and B = (l/v - 26>)/2T c 4 . Now 
we substitute these expressions into the above equation 
and integrate from a length scale Lq (temperature To) 
to a length scale L (temperature Tl) to get an equa- 
tion relating Tl to T . Then using Tl — 1/Jl gives the 
corresponding equation for Jl- Setting L/Lq = 2 n and 
29v = x, we find 



coupling throughout the low temperature phase. Now, 
for large J n and Jo, Eq. (|7|) reduces to pure power law 
behaviour viz. J 2 ~ Jq2 2 and the crossover length for 
the different temperatures can be read off as the length 
for which (Jl- Jlf 6v l{Jl -(20I/-1) J*/26v)P° v -Q be- 
comes of the same order as J 2 . Beyond the crossover 
length, we should see the correct low-temperature expo- 
nents. 
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FIG. 10. J„ (solid lines) and left hand side of Eqn.7 (sym- 
bols) as a function of system size L = 2 n for (from top to 
bottom) T = 0.3T C , 0.5T C , 0.7T C , 0.9T C . 

Thus, we have seen that the three-dimensional spin 
glass in MKA can be well described by an approxima- 
tion that is valid close to the lower critical dimension 
and that shows explicitcly that crossover scales are large. 
The analysis described above is quite general and should 
be applicable also to the Ising spin glass on a cubic lat- 
tice. There is one caveat however: the analysis we have 
carried out is valid for the couplings alone. Other quan- 
tities (for example P(q^) or -P(g' i - ) )) would have to be 
studied separately and it is possible that the crossover 
lengths would be somewhat different for different quan- 
tities. 



VI. CONCLUSION 
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where Jo is the coupling at the starting point n — and 
J n is the coupling after n iterations. 

To test whether Eq. (0) is a good approximation, we 
have inserted the values J 2 obtained from the MK recur- 
sion relations into the left hand side of Eq. ^ . We have 
chosen v = 2.8 and 9 = 0.25, in agreement with pq,p|. 
The result are the data points given in Fig. [l(]. They 
satisfy a power law with exponent 29 and a prefactor 
(Jl - J 2 ) X /(J§ - (x~l/x)J 2 )( x ~ 1 \ showing that Eq. © 
is indeed an appropriate description of the growth of the 



In this paper, we have studied the link overlap be- 
tween two identical replicas of a three-dimensional Ising 
spin glass in Migdal-Kadanoff approximation. The width 
of the link overlap distribution decreases to zero with 
increasing system size at T c as well as in the low- 
temperature phase. These findings are in agreement with 
the predictions of the droplet picture. For system sizes 
similar to the ones used in Monte Carlo simulations of a 
cubic lattice, we find the same large sample-to-sample 
fluctuations and asymmetric curve shapes as reported 
from those simulations. They must be interpreted as a 
finite-size effect and cannot be taken as an indicator for 
replica-symmetry breaking. The only reliable indicator 
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for RSB would be a width of the overlap distribution that 
does not decrease with increasing system size. 

Similarly, the link overlap in the presence of a weak 
coupling between the two replicas shows in MKA in the 
spin glass phase the singular behaviour predicted by the 
droplet picture. Data from Monte Carlo simulations are 
also in full agreement with the droplet picture. The RSB 
picture predicts a jump in the mean value of the link 
overlap at zero coupling strength that is not visible in 
Monte Carlo simulation data published so far. 

We have reproduced phenomenologically the influence 
of the critical point on the growth of the coupling con- 
stants using an approximation that is valid close to the 
lower critical dimension. This gives us a direct estimate 
of the lengthscale at any given temperature beyond which 
one needs to go in order to see zero-temperature (droplet) 
scaling without crossover effects from critical point be- 
haviour intruding too strongly. This crossover effect of- 
ten seems to get overlooked in the literature. For example 
Komori, Yoshino and Takayama |l6| in a numerical sim- 
ulation of the Ising spin glass at a temperature of 0.84T C 
found that critical scaling of the dynamics worked well: 
they found that correlation data at time t could be col- 
lapsed for systems of linear dimension L for values of L 
up to 7, by plotting against the variable t/L z ( T \ where 
z(T) is similar to the critical point dynamical exponent. 
This behaviour was interpreted by Marinari, Parisi and 
Ruiz-Lorenzo 17 as evidence against droplet scaling but 
it is clear from our work that for the system sizes studied 
at temperatures so close to T c droplet scaling would be 
quite unobservable and that critical point scaling should 
indeed work quite well. 

We therefore conclude that the droplet picture, com- 
bined with finite-size effects, can fully explain all data 
for the link overlap in the Ising spin glass. There is no 
evidence for the presence of RSB. 
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